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Identifiability of generalised Randles circuit models 

S.M.Mahdi Alavi, Adam Mahdi, Stephen J. Payne and David A. Howey 


Abstract —The Randles circuit (including a parallel 
resistor and capacitor in series with another resistor) 
and its generalised topology have widely been employed 
in electrochemical energy storage systems such as bat¬ 
teries, fuel cells and supercapacitors, also in biomedical 
engineering, for example, to model the electrode-tissue 
interface in electroencephalography and baroreceptor 
dynamics. This paper studies identifiability of gener¬ 
alised Randles circuit models, that is, whether the 
model parameters can be estimated uniquely from the 
input-output data. It is shown that generalised Ran¬ 
dles circuit models are structurally locally identifiable. 
The condition that makes the model structure globally 
identifiable is then discussed. Finally, the estimation 
accuracy is evaluated through extensive simulations. 

Index Terms —Randles circuit, Identifiability, System 
identification, Parameter estimation. 

I. Introduction 

Randles proposed an equivalent circuit for the kinetics of 
rapid electrode reactions in [T]. Since then, the model has 
been developed and has become the basis for the study 
of many electrochemical energy storage systems such as 
batteries, fuel cells and supercapacitors, [2], [3], Figure □ 
shows a generalised Randles model consisting of an ohmic 
resistor, Roc, in series with a number of parallel resistors 
and capacitors, and a capacitor C w . In electrochemical 
applications, the ohmic resistor represents usually 
the conduction of charge carriers through electrolyte and 
metallic conductors. The resistors and capacitors in the 
parallel pairs represent the charge transfer resistance and 
the double layer capacitance, respectively, or are an ap¬ 
proximation of a diffusion process, a. 0 - The number 
of parallel R 's and C s depends on how many of these 
pairs are required such that the frequency response of the 
generalised Randles model fits with the device impedance 
spectra within the frequency range of interests 0-0- 
For instance, in [3], the number of the parallel pairs was 
determined by minimising the error between the model 
and measured voltages. The capacitor C w , also known as 
the Warburg term, accounts for a diffusion process, HI- 0; 
or in a battery or supercapacitor, it may represent state of 
charge [9]. It should be noted that the open-circuit voltage 
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source is not considered in the generalised Randles model 
Figure jT] and it is only focused on the impedance model 
in this paper. 

The generalised Randles circuit models have also been 
employed in biomedical engineering. An electric cir¬ 
cuit model of the electrode-tissue interface in electroen¬ 
cephalography includes two parallel R-C pairs in series 
with two resistors and two voltage sources m, which is a 
special case of the circuit given in Figure[Q It is noted that 
C w can be considered as a parallel R and C with R = oo. 
The viscoelastic analog of the generalised Randles model is 
employed in cardiovascular and cerebral haemodynamics 
modeling, in describing the viscoelastic properties of the 
aortic wall and the coupling of the nerve endings of the 
baroreceptor neurons in the carotid sinus and aortic arch 
nu. m and in relating the fluctuation of the arterial 
blood pressure with the cerebral blood flow velocity ESI- 

The identification or parameter estimation of the gen¬ 
eralised Randles model (with different numbers of parallel 
R’s and C’s) is important in condition monitoring, fault 
diagnosis and control [3], 0, [13]— (2Dj. In [T3] , the authors 
showed that the Randles circuit can be used for monitoring 
the battery charge transfer overvoltage. In Ha. identifica¬ 
tion tests for parameter estimation of lead-acid batteries 
are suggested. In [2D], identification through fitting the 
impedance spectra in the frequency domain is presented. 

The objective of this paper is to study the identifiability 
of the generalised Randles model shown in Figure [IJ 
that is, whether the model parameters can be estimated 
uniquely from input-output data. Typically the identifia¬ 
bility problem is divided into two broad areas: parameter 
estimation accuracy and structural identifiability. Param¬ 
eter estimation accuracy considers the practical aspects 
of the problem that come with real data such as noise 
and bias 1211 . In studying structural identifiability, on the 
other hand, one assumes that noise-free informative data 
is available and therefore it is, in fact, a data-independent 
concept. Unidentifiable parameters can be assigned an 
infinite number of values yet still lead to identical input- 
output data. Thus, structural identifiability is a necessary 
condition for identifiability and parameter estimation. A 
number of analytical approaches to structural identifia¬ 
bility have been proposed, including Laplace transform 
(transfer function) [23, [23], Taylor series expansion [211 . 
[23] . similarity transformations [2Bj-[33], and differential 
algebra [31], [33] ■ I n linear systems, it has been shown 
that controllability and observability properties are closely 
related to the concept of structural identifiability, 03- 
m- For example, it was shown that single-input single¬ 
output linear time-invariant systems are structurally iden¬ 
tifiable if and only if their observer canonical form is 


2 


Ra 


Rl R 2 Rn 

— / WW'— —'WW'— r^MA^-i 


^WWH 




<7 


C 2 


C n 


C w i(t) 


+v w 


+V!~ 


+V 2 - 


+v n - 


Fig. 1. The generalised Randles equivalent circuit model. 


controllable (see Chapter 4 in [33]). However, controllable 
and observable systems can still be unidentifiable in the 
general case m- 

Recently, there has been a significant interest in the 
identifiability analysis of battery models [SB -[33]. In m, 
the structural identifiability of a five-element equivalent 
circuit model (ECM) including two capacitors was dis¬ 
cussed by comparing the number of unknown parameters 
of the Transfer Function (TF) and the circuit. In [37], the 
structural identifiability of a more general nonlinear ECM 
is analysed based on the observability conditions. It is 
shown that cells with serial connections are not observable, 
demonstrating that complete estimation of the state-of- 
charge and model parameters from lumped measurements 
with series cells is not possible trough independent mea¬ 
surements. However, it is shown that lumped models with 
parallel connectivities are observable provided that none 
of the parallel cells are identical. The identifiability of 
battery electrochemical models was discussed in m-m- 
In particular, it was shown that some of the electrochem¬ 
ical model parameters are not identifiable given typical 
charge-discharge cycles [3H]-[1D]- In [H], it was shown that 
the shape of the charge-discharge cycles plays a crucial 
role in the identifiability of battery parameters. It was 
also demonstrated that the system identification method 
can be employed in monitoring the battery film growth. 
The identifiability problem of Randles ECMs has been 
studied in [32] . [33], where the authors consider models 
that include up to two capacitors, and the analysis is 
based on the Fisher information matrix (FIM). The FIM 
provides some information about the sensitivity of the 
measurement to the model parameters by using likelihood 
functions, [44] . In [42], a bound of estimation errors was 
developed by using the Cramer-Rao theorem. In [35] . 
a method was proposed to optimally shape the battery 
cycles and improve the identifiability. 

This paper shows that the generalised Randles model 
in Figure [T| is structurally globally identifiable for n = 
1 , and structurally locally identifiable for any finite n > 
1 , and becomes globally identifiable assuming an ordering 
through the generalised circuit. Finally, the identifiability 
of the model is assessed through extensive simulations. 

II. The model parameterisation and problem 

STATEMENT 

The state-space and TF parameterised models of the 
generalised Randles circuit given in Figure [T] are derived. 


Define the electric current as the system input u(t) = 
i{t) £ R, the terminal voltage as the system output y(t) = 
v(t) £ R, the voltages across the internal R-C pairs as the 
states x(t) = [ vi(t) ■■■ v n (t) v w (t) ] T £ R” +1 , and 
the model parameters as 

0=[R oo Rl ■■■ R n Cl ■■■ Cn C w ] , (1) 

where 6 belongs to some open subset T> C R 2n+2 . 


A. State-space parametrisation 

By using Kirchhoff’s laws, a state-space model structure 
of the Randles circuit, parameterised by 6 , can be written 
as: 

f ^ x 0) = MO) x(t) + B{6) u(t ) ^ 

1 y(t) = C{0) x(i) + D{9)u(t) 
where A(0) £ R(«+i)x(n+i) ) e Rn+i, c ^t g R n+ 1 

and D(0) £ R are matrices that depend on the parameter 
vector 6 , and are given by 
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, D{<S) = d, (3) 


with time constants Ti = RiCi for i = 1, • • • , n and 

5 bi — 77 > b w — — . d — Rco • (4) 

7*7, C^7/; 


B. Transfer function parametrisation 


Denote the model’s TF by T(s,6), where s represents 
the Laplace operator. Using the formula T{s,0) = C{sl — 
A)~ 1 B-\-D, a parameterised TF of the generalised Randles 
circuit can be written as: 


T{s,0) 


bi{0) b w {0 ) 

“ s + ai{0) s 
1=1 


+ D(6). 


(5) 
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C. Problem statement 

Determine the conditions for the model structure ED 
(equivalently in the state-space form ED and ED): where 8 
is an unknown parameter vector ED, to be locally and/or 
globally structurally identifiable. 


III. Main results 

Following [33] a definition of structural identifiability is 
given as follows. 

Definition 1: Let Ad be a model structure with the TF 
T(s,8), parametrized by 8 , where 8 belongs to an open 
subset T>t C R m , and consider the equation 


T(s,0) = T(s,8*) for almost all s, (6) 


where 8,8* £ D T . Then, the model structure Ad is said to 
be 


- globally identifiable if El has a unique solution in T>t, 

- locally identifiable if ED has a finite number of solu¬ 
tions in T>t, 

- unidentifiabile if ED has a infinite number of solutions 
in T>t- 

Remark 1: Instead of Definition [L] one can use the 
so-called coefficient map defined as follows, [25] . 
Consider the monicQ TF: 


T(s,8) = 


Co(9) + Ci(6)s ■ 


+ c fei ( 8)s kl 


do(8) + di(0)s H-b d k2 -i{8)s k2 ~ l + s k2 ' 


(7) 

and associate with it the following coefficient map Ct : 
R m D T> t —t R fel+fc2+1 defined as 


Ct{ 8) — c 0 {6),. . ■, Cfcj (8),do{0), • • -, dk 2 -i{8) 


( 8 ) 


The model structure Ad is globally identifiable if the coef¬ 
ficient map Ct is one-to-one (injective); locally identifiable 
if Ct is many-to-one; and unidentifiable if it is infinitely 
many-to-one. 

The following lemma will be used in the proof of the 
identifiability of the generalised Randles circuit model. 

Lemma 1: Let (7t(to) be a coefficient map associated 
with the TF 


T(s,8) 


E 


bi 


s + a* 


+ d, 


where 8 — (ai,..., a m , b ±,..., b m , d) and a, £ R for 
i = 1,..., to are pairwise different. Then the following 
statements hold: 

(a) If to = 1, then Ct (to) is one-to-one. 

(b) If to > 1, then Ct [in) is many-to-one. 

Proof. Part(a). For to = 1, the identifiability equation 
ED is given by 


s + Oi S-(-ftj 

and has a unique solution (ai,&i,d) = (aj,6|,d*), which 
proves part (a). 


1 In a monic TF, all coefficients are normalised such that the 
coefficient of the greatest order in the denominator is 1. 


Part(b). For to > 1, the identifiability equation ED is 
given by 


E 


bj 

s + ai 


+ d = 


Z^ 


+ d*. 


(9) 


We claim that equation ED admits only finite (more 
precisely, to! = 1 x 2 x • • • x to) number of solutions. To 
prove the claim note that ED is the equality of two rational 
functions, which is satisfied provided that 


(s + ai) x • • • x (s + a m ) = (s + a*) x ■ • • x (s + a* m ). (10) 


Since m distinct roots uniquely characterise a monic 
polynomial of degree to, and there are to! permutations 
of n roots of (s + ai) x ••• x (s + a m ), equation m 
has to! solutions. Now, let us fix the permutation s : 
{l,...,m} —> {1 ,...,to} and consider an assignment 
(ai,... ,a m ) = (a* ,..., a* ). Since a* for i = 1 ,...,m are 
assumed to be pairwise distinct, the expressions 1 /(s+aj), 
thought of as functions of the variable s, are linearly 
independent. Finally, since each side of equation ED is a 
linear combination of linearly independent functions, we 
immediately obtain that bj = b* for i = 1 ,...,m. This 
concludes the proof of part (b). □ 

Now, we introduce the concept of reparametrisation (see 
e.g. [28] h which we will use the proof of our main theorem. 


Definition 2: A reparametrisation of the model structure 
Ai with the coefficient map Ct is a map 1Z : R fc D Dt —>■ 
R m such that 

Im (Ct 0 7Z) = IvhCt , (11) 

where Im denotes the image of the map. Moreover, the 
reparameterisation is identifiable if the map Ct ° TZ '■ R fc A 
D —> R fc i+fc 2 +i j s identifiable. 

The main result of this section is the following theo¬ 
rem which describes the identifiability of the generalised 
Randles circuit. 

Theorem 1: Let M. rc(ji) denotes the state-space model 
structure ED with the matrices ED parametrised by ED, 
where n is the number of parallel RC elements connected 
in series (see Figure [1}. Then the following conditions hold: 

(a) If n = 1, then the model structure AiRc(n) is globally 
identifiable. 

(b) If n > 1, then the model structure AiRc(n) is locally 
identifiable. 

(c) If n > 1, and there is an ordering through the 
generalised circuit as 

a n ^ Cn—i ^ ^ iR, (12) 

then the model structure Ad Rc(n) is globally identi¬ 
fiable. 

Proof. Consider the model structure Ai rc(u) and let 
T(s,8) denote the corresponding TF given by ED: where 
the parameter vector 8 is given by (Q. We write T(s,8) 
as a rational function ED of degree Aq = k% = n + 1 and 
associate with it the coefficient map C' T : R 2n+2 D V — > 
R 2n + 3 . 
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Part (a). For n = 1, the coefficient map C' T : R 4 D T> —> 
R 5 can be written explicitly as 


C T (0) 


c o (0),c 1 (O),c 2 (O),d o (0),d 1 (O) , 


where 

com = 


i 


RiC\C u 


ci{0) = 


RooCw -fi R\C W + R\C\ 


c 2 {0) = Roo, do(0) = 0 , d i( 0 ) = 


/i’i C Y | c u 

1 


R\Ci 


By direct computation, we can check that equation 
C T {0) = C T (0*) admits a unique solution 0 = 0*. Thus 
the coefficient map is one-to-one, and the model structure 
is globally identifiable. 

Part (b). For n > 1, the coefficient map C T can be 
written as the following composition 


C t =C t o TZ c , 


(13) 


where the map 1Z C : R 2n+2 D V —> R 2n+2 is the 
reparametrisation 1Z C (0) = (ai,..., a n , b±,..., b n , b Wl c?) 
defined by 0; and Ct is the coefficient map associated 
with the TF 

n , , 

T(s, 0 a b) = —y-1—— + d, (14) 

^ s + ai s 
1— 1 


where 0 afi = (ai,. •., a n , bi, ..., b n , b w , d) G R 2n+2 . By 
Lemma [Tj (for m = n + 1 and a n+ 1 = 0) the map Ct is 
many-to-one; and 1Z C is a one-to-one map with an inverse 


Roo — d , C w — , Ri , Ci — , i — 1 ,..., ti. 

Uyj Q-i Oi 


The map C T is many-to-one, since it is a composition of 
a one-to-one with many-to-one map. Thus for n > 1, the 
model structure M.R.c{n) is locally identifiable. 

Part (c). Finally, the identifiability equation © under 
the condition m admits a unique solution, which con¬ 
cludes part (c). □ 

Remark 2: The same procedure is applicable to the 
discrete time model. The Euler’s first order approximation 
is the simplest approximation. Its identifiability analysis 
is easier because the coefficients of the discrete time TF 
equals the number of parameters. However, it might lead 
to numerical instability. If higher order approximations are 
applied, the challenge remains to prove whether there is 
a finite-to-one maps between the discrete time model and 
parameters. 


IV. Simulations 

Consider a 6 -element Randles circuit with R^, C w Ri , 
Q, i = 1,2. In this section, accuracy of the estimation 
in the presence of noise-free and noisy data, subject to 
random initial guess of estimations is studied. First the 
informative data set used in the simulations and the way 
it is generated are described. More details are provided in 
references m, m-m- 


A. Generation of informative data set 


A data set is informative if the input is persistently 
exciting. A persistently exciting input adequately excites 
all modes of the system. In linear systems, the order of 
the system determines the order of persistent excitation. 
The order of a persistent excitation equals the number of 
coefficients of the monic TF that need to be identified, (see 
Theorem 13.1. in [32]). The monic TF of the 6 -element 
Randles circuit is 


/ 3 s 3 + f 2 s 2 + fis + fp 
s 3 + g 2 s 2 + gis + g Q 


(15) 


with unknown coefficients to be identified. Therefore, the 
necessary order of the persistently exciting input is 7. In 
the frequency domain, this means that the spectrum of 
the excitation input should have at least 7 nonzero points. 
This simulation focuses on multi-sine excitation signals 
which are widely employed in electrochemical impedance 
spectroscopy (EIS) techniques, [3], [1], [IS]. The same 
procedure can be applied to other inputs. The multi-sine 
signal is given by 


1 

u(t) = y, nij cos (ujt + ifj), (16) 

3 =1 


where l represents the number of sinusoids, and mj, ujj 
and <pj G [—7r, 7 r) denote the magnitude, frequency in 
radians per second and phase in radians, respectively. The 
spectrum of the multi-sine signal is given by 


ffl- 

$ u (w) = 2 tt^-A[<5(w - Wj)+ <5(w + Wj)], (17) 

3= 1 

where S(cu) is the delta function or impulse at frequency 
ujj. The spectrum of each sinusoid signal contains two 
nonzero points, therefore 4 sinusoid signals are enough 
to generate the informative data set for the 6 -element 
Randles circuit model. 


The magnitudes mj, frequencies ojj and phases 4>j are 
arbitrary real values. Specific applications might impose 
additional constraints. For instance, in EIS techniques, the 
magnitude of the input signal may vary from milliampere 
to ampere depending on the size of the energy storage 
system. For the sake of simplicity, all magnitudes are 
assumed equal to mj = 10 -3 for j = 1, ,4. The 

frequencies could be equally or logarithmically spread over 
the frequency band w m i n < w < w max . The values of w m i„ 
and w max depend on the system dynamics. If Wj’s are 
equal, the Schroeder phase choice is suggested to reduce 
the Crest factor 33]. 1 43] , 

Remark 3: A multi-sine signal may not exactly be zero 
mean. Sometimes a uni-directional current flow is required. 
In this case, the multi-sine excitation current needs to be 
superimposed on a known constant DC offset current, (see 
examples in [3] and [43] ). In system identification, the data 
are typically pre-filtered to remove these types of offsets. 









5 



(a) Excitation input u(t). 


Fig. 2. The multi-sine input excitation signal a nd its associated output 
and C w = 300 F. The voltage response, Figure|2(b)| is computed using 
a = 10“ 4 . The DC offset has been removed to have zero-mean signals 



(b) Associated output y(t). 


for Roo = 0.05 $2, Ri = 0.2 ft, Ci = 0.3 F, R 2 = 0.4 ft, C 2 = 0.6 F, 
the model © in the noise-free and noisy case with standard deviation 


B. The results discussion 


The following operating point is arbitrarily selected [3], 

01 : 


Roo = 0.05 n, Rx = 0.2 Cl, Cl = 0.3 F 
i?2 = 0.4 Cl, C 2 = 0.6 F, C w = 300 F 


The multi-sine excitation signal, Figure 2(a) is generated 
by using model (flfil) with l = 4, nij = 10 -3 , f m i n = 0.2 
Hz, / m ax = 500 Hz and Schroeder phase with </i randomly 
chosen at 1.9775. The Crest factor of the signal is 1.9999. 
The voltage response, Figure |2(b)[ has been computed 
using the model ©. The DC offset has been removed to 
ensure zero-mean signals. 


Given the true values, the smallest time constant is 
T m i n = i?iCi = 0.06 s. The sampling time should be 
several times larger than l/r min . In this simulation, the 
sampling frequency is chosen at f s = 500 Hz, which is 
16.6 times greater than the inverse of the minimum time 
constant. The test duration should typically be several 
times larger than the maximum time constant, which is 
Tmax = RooC w = 15 s. In 07 ] , 6 to 8 times r max is 
suggested, however, this might vary in different applica¬ 
tions. In this simulation, a test duration of td = 100 s is 
applied. The continuous-time TF (fl 5 l) is identified using 
Matlab’s system identification toolbox m- The circuit 
parameters are calculated directly from the coefficients of 
the TF using the formulas shown in Table II and discussed 
later. In order to study the consistency of results, each 
test is repeated 100 times, every run with a random initial 
guess of parameters. The roots of the denominator must 
be positive real numbers. Those estimations, which lead to 
complex or negative poles for the TF, or to C w > 1000 or 
to Ci > 10, i = 1,2, are considered outliers and discarded 
from the analysis. 


The relative mean error, e r is defined as follows: 

true value — mean of estimations 


e r = 100 x 


true value 


%. (18) 


The estimation accuracies from noise-free and noisy 
data with a zero-mean noise and standard deviation a = 
10~ 4 are compared together. The signal to noise ratio is 
= 10. Table U shows the mean, standard deviation 
and relative errors of the estimations. Figure [3] shows the 
histograms of the estimations from the noisy data. 

Regardless of the random initial condition, in more than 
80% of simulations, the noise-free estimations converge 
with e r less than 10% for R i: Ci, i = 1,2 and for C w . 
The largest e r is 10.36% for R^ (Table HJ). 

Figure [3] shows that estimations from noisy data are 
more distributed around the true values. However as it is 
seen in Table U the mean values of estimations in both 
noise-free and noisy cases remain the same. The relative 
mean errors of estimations are again less than 10% for all 
parameters except of R^, which is again the largest and 
10.38%. Also, e r ’s increase with noisy data compared to 
that of in the noise-free case. 

The largest standard deviation is for C w . This could be 
because of the pure integrator associated with C w , which 
appears in the TF, which might require some modifications 
of the data set, m ■ For instance, a methodology has 
been proposed in .36] that removes the integral term by 
modifying the input signal as J m odified(f) = f 


C. Calculating the coefficients of R — R\\C — R\\C — C 
circuit 

Here we show how to compute the parameters of the 
6-element Randles circuit. Using (0, the circuit’s TF is: 

u, /3s 3 + f 2 s 2 + f\S + /o 
s i + g 2 s z + gis + g 0 
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Fig. 3. Histogram of accepted estimations in 100 runs from noisy 
data. 


The relationships between the coefficients and the circuit 
parameters are given by: 

/3 = d, / 2 = b\ + b-2 + b w + (ai + a,2)d 

fi = a 2 &i + aib ’2 + (ai + a 2 )b w + fo = a\a 2 b w (19) 

gi = ai + 0,2, gi = aia2, 3o = 0 

where the parameters ai, bi, b w and d are as defined in 

©• 

Because the circuit has an integrator, the identification 
method should be set up such that a pole of the denomina¬ 
tor is fixed at s = 0. The identification software typically 
allows to fix a number of poles and zeros at certain values, 
m- Using the first equation of ITT71) and (J3D , Roo is given 
by: 

Roo=h 

The roots of s 3 + g 2 S 2 + g\S + go = 0 are a\ and 0 , 2 , 
and the one that has been fixed at s = 0. By using the 
condition m, select the smallest root as 0 , 2 , and the 
remaining root as ai. From (JTj), the circuit’s time constants 
are obtained as follows: 

Ti — — for i — 1,2. 
ai 

From (fl9l) and (fT^l) ,b\,b 2 and b w are obtained by solving 


the following set of equations for X: 


1 

1 

1 


f 2 - (ai + a 2 )/3 

a 2 

a\ 

ai + a 2 

X = 

fi ~ aia 2 / 3 

0 

0 

a\ a 2 


fo 


where, 

X = [ bi &2 b w ] 


Other parameters of the circuit are subsequently ob¬ 
tained as: 


C u 


Ci = 


b<’ 



i = 1,2. 


The parameters of different topologies can simply be 
calculated using the same approach. Table El provides 
formulas for four widely used Randles models. 


V. Conclusions 

We showed that the generalised Randles circuit model 
is locally identifiable and the model structure becomes 
globally identifiable if an ordering through the circuit is 
assumed. The results were confirmed through extensive 
simulations. Finally, explicit formulas for the coefficients 
of widely used Randles circuits were presented. 
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TABLE I 

Estimation results. 



Parameter 

Roo 

Ri 

Ci 

R-i 

c 2 

c w 


True Value 

0.05 

0.2 

0.3 

0.4 

0.6 

300 


mean 

0.0552 

0.1883 

0.2919 

0.4061 

0.5866 

289.3709 

Noise free (9 outlier estimates) 

st.d. 

0.0010 

0.0393 

0.0617 

0.0360 

0.0734 

48.3550 


e r 

10.36 

5.87 

2.68 

1.52 

2.22 

3.54 


mean 

0.0552 

0.1847 

0.2886 

0.4094 

0.5803 

300.9287 

Noise <7 = 1(T 4 (11 outlier estimates) 

st.d. 

0.0012 

0.0451 

0.0706 

0.0415 

0.0867 

87.8931 


e r 

10.38 

7.63 

3.79 

2.34 

3.28 

0.31 


TABLE II 

The coefficients of four widely used Randles circuit models. 










































